Parasitic small-moment-antiferromagnetism and non-linear coupling of hidden order 
and antiferromagnetism in URu2Si2 observed by Larmor diffraction 
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We report simultaneous measurements of the distribution of lattice constants and the antiferro- 
magnetic moment in high-purity URu2Si2, using both Larmor and conventional neutron diffraction, 
as a function of temperature and pressure up to 18 kbar. We establish that the tiny moment in the 
hidden order (HO) state is purely parasitic and quantitatively originates from the distribution of 
lattice constants. Moreover, the HO and large-moment antiferromagnetism (LMAF) at high pres- 
sure are separated by a line of first-order phase transitions, which ends in a bicritical point. Thus 
the HO and LMAF are coupled non-linearly and must have different symmetry, as expected of the 
HO being, e.g., incommensurate orbital currents, helicity order, or multipolar order. 
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In recent years hydrostatic pressure has become widely 
used in the search for new forms of electronic order, 
because it is believed to represent a controlled and 
clean tuning technique. Novel states discovered in high- 
pressure studies include superconducting phases at the 
border of magnetism and candidates for genuine non- 
Fermi liquid metallic states. However, a major un- 
certainty in these studies concerns the possible role of 
pressure inhomogeneities, that originate, for instance, in 
the pressure-transmitting medium and inhomogeneities 
of the samples. To settle this issue requires microscopic 
measurements of the distribution of lattice constants 
across the entire sample volume, which, to the best of 
our knowledge, has not been available so far. 

The perhaps most prominent and controversial exam- 
ple that highlights the importance of sample inhomo- 
geneities and pressure tuning is the heavy-fermion su- 
perconductor URu2Si2. The reduction of entropy at a 
phase transition at Tq ~ 17.5 K, discovered in URu2Si2 
over twenty years ago, is still not explained The 
associated state in turn is known as 'hidden order' (HO). 
The discovery of the HO was soon followed by the ob- 
servation of a small antiferromagnetic moment (SMAF), 
nis ~ 0.01 — 0.04 /iB per U atom [4] then believed to be 
an intrinsic property of the HO. The emergence of large- 
moment antiferromagnetism (LMAF) of ~ 0.4 /iB 
per U atom [5] under pressure consequently prompted 
intense theoretical efforts to connect the LMAF with the 
SMAF and the HO. In particular, models have been pro- 
posed that are based on competing order parameters of 
the same symmetry, i.e., linearly coupled order parame- 
ters, in which the SMAF is intrinsic to the HO 0,0, 11,0 • 
This is contrasted by proposals for the HO parameter 



such as incommensurate orbital currents jlOl], multipolar 
order [ll[, or helicity order pJi], where HO and LMAF 
break different symmetries. 

The symmetry relationship of HO and LMAF clearly 
yields the key to unravelling the nature of the HO state 
[1, 0]. While some neutron scattering studies of tfie 
temperature-pressure pfiase diagram suggest tfiat tfie 
HO-LMAF phase boundary ends in a critical end point 
, other studies concluded that it meets the boundaries 
of HO and LMAF in a bicritical point [13, E, El, E^- 
This distinction is crucial, as a critical end point (bicriti- 
cal point) implies that HO and LMAF have the same (dif- 
ferent) symmetries, respectively [9]. Moreover, there is a 
substantial disagreement w.r.t. the location and shape of 
the HO-LMAF phase boundary (see e.g. Ref. [I^). This 
lack of consistency is, finally, accompanied by consider- 
able variations of the size and pressure dependence of the 
moment reported for the SMAF ^] , where NMR studies 
in p owder samples suggested the SMAF to be parasitic 
[l9|. Accordingly, to identify the HO in URu2Si2 it is es- 
sential to clarify unambiguously the nature of the SMAF 
and the symmetry relationship of HO and LMAF. 

It was long suspected that the conflicting results are 
due to a distribution of lattice distortions due to de- 
fects. Notably, uniaxial stress studies showed that LMAF 
is stabihzed if the c/a ratio r] of the tetragonal crys- 
tal is reduced by the small amount Ar]c/r] ~ 5 • 10~^ 
(2q| . Hence, the parasitic SMAF may in principle re- 
sult from a distribution of r] values across the sample, 
with its magnitude depending on sample quality and ex- 
perimental conditions. In particular, differences of com- 
pressibility of wires, samples supports or strain gauges 
that are welded, glued or soldered to the samples will 
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total Larmor phase of precession ^ thereby depends hn- 
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FIG. 1: (a) Schematic of Larmor diffraction [22. [23l|. see text 
for details, (b) Typical variation of the polarization P as a 
function of the total Larmor phase (c) Pressure depen- 
dence of the width of the lattice- const ant distribution for the 
a- and c-axis in URu2Si2. With increasing pressure the width 
of the distribution increases. 



forcibly generate uncontrolled local strains that strongly 
affect any conclusions about the SMAF signal (see, e.g., 

Refs [13, EH [13) • 

In this Letter we report simultaneous measurements of 
the distribution and temperature dependence of the lat- 
tice constants, as well as the antiferromagnetic moment, 
of a pure single crystal of URu2Si2 utilizing a novel neu- 
tron scattering technique called Larmor diffraction (LD). 
This allowed to study samples that are completely free to 
float in the pressure transmitting medium, thereby expe- 
riencing essentially ideal hydrostatic conditions at high 
pressures. Our data of the distribution of lattices con- 
stants f{Ar]/r]) establishes quantitatively that the SMAF 
is purely parasitic. In addition, we find a rather abrupt 
transition from HO to LMAF which extends from T = 
up to a bicritical point (preliminary data of Tn{p) were 
reported in [lH). We conclude that the HO-LMAF tran- 
sition is of first order and that HO and LMAF must be 
coupled non-linearly. This settles the perhaps most im- 
portant, long-standing experimental issue on the route to 
identifying the HO. 

Larmor diffraction permits high-intensity measure- 
ments of lattice constants with an unprecedented high 



resolution of Aa/a ^ 10"^ [22|,l23|. As shown in Fig.[T](a) 



the sample is thereby illuminated by a polarized neutron 
beam (arrows indicate the polarization); G is the recipro- 
cal lattice vector; 6^^ is the Bragg angle; AD is the polar- 
ization analyzer and detector. The radio frequency (RF) 
spin resonance coils (C1-C4) change the polarization di- 
rection, as if the neutrons undergo a Larmor precession 
with frequency col along the distance L = Li + L2. The 



the mass of the neutron [22|, |23|). 



The Larmor diffraction was carried out at the spec- 
trometer TRISP at the neutron source FRM II. The tem- 
perature and pressure dependence of the lattice constants 
was inferred from the (400) Bragg peak for the a axis and 
the (008) Bragg peak for the c axis. The magnetic or- 
dered moment was monitored with the same setup using 
conventional diffraction. For our high-pressure studies a 
Cu:Be clamp cell was used with a Fluorinert mixture [23 ]. 
The pressure was inferred at low temperatures from the 
(002) reflection of graphite as well as absolute changes 
of the lattice constants of URu2Si2 taking into account 
published values of the compressibility. 

The single crystal studied was grown by means of an 
optical floating-zone technique at the Amsterdam/Leiden 
Center. High sample quality was confirmed via X-ray 
diffraction and detailed electron probe microanalysis. 
Samples cut-off from the ingot showed good resistance 
ratios (20 for the c axis and ~ 10 for the a axis) and a 
high superconducting transition temperature Tc ~ 1.5 K. 
The magnetization of the large single crystal agreed very 
well with data shown in Ref. [25[ and confirmed the ab- 
sence of ferromagnetic inclusions. Most importantly, in 
our neutron scattering measurements we found an anti- 
ferromagnetic moment ~ 0.012 /hb per U atom, which 
matches the smallest moment reported so far |l8| . 

As the Larmor phase <l> is proportional to the lattice 
constant a, the polarization P of the scattered neutron 
beam reflects the distribution of lattice constants across 
the entire sample volume [2^. While Larmor diffraction 
was recently employed for the first time in a high-pressure 
study [26], measurements of the distribution of lattice 
have not been exploited to resolve a major scientific issue 
(for proof of principle studies in Al-alloys see Ref. [22|). 

In order to explore the origin of the SMAF moment we 
have measured the spread of lattice constants, keeping in 
mind that a distribution of the c/a lattice-constant ratio 
T] may be responsible for AF order in parts of the sam- 
ple. As shown in Fig.[T](b) P(^) changes only weakly be- 
tween ambient pressure and 17.3 kbar. Accordingly the 
distribution of lattice constants, which is given by the 
Fourier transform of P{^)^ changes only weakly as a func- 
tion pressure as shown in Fig. [lb). Thus pressure only 
slightly boosts the distribution of lattice constants, but 
does not generate substantial additional inhomogeneities. 
Moreover, we confirmed that the size of remained un- 
changed tiny after our high-pressure studies. 

Assuming a Gaussian distribution of both lattice con- 
stants, we infer the distribution of 77. At low pres- 
sure, we arrive at a full width at half-maximum of 
f{Ar]/r])YwuM ^ 3.8 • 10~^. Recalling that the tail of the 
Gaussian distribution beyond Ar]c/v ^ 5 -10"^ represents 
the sample's volume fraction in which LMAF forms [2q| . 
an average magnetic moment of 0.015(5)/iB is expected in 
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FIG. 2: Temperature dependence of the lattice constants and 
thermal expansion at various pressures. The HO and LMAF 
transitions are indicated by empty and filled arrows, respec- 
tively, (a) Data for the a-axis; (b) thermal expansion of the 
a-axis derived from the data shown in (a), (c) Data for the 
c-axis; (d) thermal expansion of the c-axis derived from the 
data shown in (c). 
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FIG. 3: (a) Pressure dependence of the low-temperature mag- 
netic moment. The inset shows the temperature dependence 
of the magnetic peak height at various pressures, (b) Phase 
diagram based on Larmor diffraction and conventional mag- 
netic diffraction data. The onset of LMAF is marked by full 
and of HO by empty symbols (x marks a transition near base 
temperature) . For better comparison data of Tn (black sym- 
bols) from Refs. [3, El, 113, El are shown. The To values (red 
symbols) of all references are consistent. Inset: background- 
free specific-heat jump derived from thermal-expansion data 
via the Ehrenfest relation (full circles) . Heat capacity data is 
taken from Refs. [1| (square) and [3 (empty circle). 



our sample. This is in excellent quantitative agreement 
with the experimental value and represents the first main 
result we report. 

We continue with a discussion of the temperature de- 
pendence of the lattice constants and its change with 
pressure. At zero pressure in a wide temperature range 
below room temperature (not shown) the thermal expan- 
sion is positive for both axes [l^. However, for the c-axis, 
the lattice constant shows a minimum close to 40 K and 
turns negative at lower temperatures. This general be- 
havior is unchanged under pressure. 

Shown in Fig. [2] are typical low-temperature data for 
the a- and c-axis. At ambient pressure Tq can be barely 
resolved. However, when crossing pc ^ 4.5 kbar a pro- 
nounced additional signature emerges rapidly and merges 
with To . Measurements of the antiferromagnetic moment 
(see below), identify this anomaly as the onset of the 
LMAF order below its T^. At the transition the lattice 
constant for the a- and c-axis show a pronounced con- 
traction and expansion, respectively. 

The pressure and temperature dependence of deter- 
mined at (100) is shown in the inset of Fig.[3](a). Close to 
Pc ^ 4.5 kbar, where Tn is near base temperature, we ob- 
serve first evidence of the LMAF magnetic signal, which 



rises steeply and already reaches almost its high-pressure 
limit at 5 kbar. The pressure dependence was determined 
by assuming the widely reported high-pressure value of 
rris = 0.4/iB and comparing the magnetic (100) and nu- 
clear (004) peak intensity. 

The temperature-pressure phase diagram shown in 
Fig.[3](b) displays Tq and Tat taken from the magnetic 
Bragg peak (Fig.[3](a)) and from the Larmor diffraction 
data (Fig. [2]), respectively. The different data sets show 
excellent agreement (except at p = 0, due to the variable 
parasitic nature of the SMAF). The main results con- 
tained in Fig. [3] are (i) the very small value of 0.012 fiB 
of the average low-temperature ordered moment at zero 
pressure, (ii) a particularly abrupt increase of the low- 
temperature moment at Pc ^ 4.5 kbar (as compared to 
previous studies BEE EH), (iii) the steep slope 
of the HO-LMAF phase boundary, and (iv) the merging 
of the HO-LMAF phase boundary with the Tq transition 
lines at approximately 9 kbar. Fig. [2](b) shows that we 
can follow this phase boundary from low T up to Tq. 

Most importantly, (iv) implies that HO and LMAF are 
fully separated by a phase boundary which has to be of 
first order with a bicritical point, since three second-order 
phase transition lines cannot meet in one point. (Note 
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FIG. 4: Extended phase diagram of URu2Si2 based on the 
data presented here (diamonds) and signatures in resistivity 
(brown[2^ and grey[3 circles). Tmax o and Tmax,Mc denote 
coherence maxima in the resistivity \2§\ and magnetisation 
[25], respectively. The HO might either mask a QCP (I) or 
replace LMAF (II) near quantum criticality. 



that the phase boundaries from the HO and LMAF to the 
disordered high-temperature phase are already known to 
be of second order from qualitative heat capacity mea- 
surements [ll,[l6|.) This conclusion is perfectly consistent 
with (ii) and (iii), and it excludes a linear coupling be- 
tween the HO and LMAF order parameters [9|. 

Although our results show that the HO and LMAF 
must have different symmetry, they also suggest that 
both types of order may have a common origin. Using 
the Ehrenfest relation, we have converted our thermal 
expansion data into a background-free estimate of the 
specific-heat jump, AC, at the PM to LMAF transition. 
The inset in Fig.[3](b) shows a continuous evolution with 
pressure of the jumps at the LMAF to PM and corre- 
sponding HO to PM transitions. The common origin of 
both phases may in fact be related to excitations seen in 
neutron scattering. Notably, excitations at (1,0,0) only 
appear in the HO phase and may be its salient prop- 



erty 2^, |30|. However, excitations at (1.4,0,0) in the PM 
state become gapped in both the HO and LMAF state 
and have been quantitatively linked to the specific heat 
jump at zero pressure [sij . 

Finally, Fig. [H suggests a new route to the HO. It shows 
a summary of the pressure evolution of features in the re- 
sistivity and magnetization, which suggest the existence 
of an AF quantum critical point (QCP) at an extrap- 
olated negative pressure. Remarkably, these features, 
which are usually denoted as Kondo or coherence temper- 
ature, all seem to extrapolate to the same critical pres- 
sure. The HO in URu2Si2 could be understood as emerg- 
ing from quantum criticality, possibly even masking an 
AF QCP (case I in Fig.[4j). However, considering the 
remoteness of the proposed QCP, the balance between 
LMAF and HO might rather be tipped by pressure in- 
duced changes of the properties of URu2Si2. This has 
for example been suggested in a recent proposal, where 



HO and LMAF are believed to be variants of the same 
underlying complex order parameter ^2[. In such a sce- 
nario, the HO is more likely to collapse at the proposed 
QCP as denoted by the dashed line (case II in Fig.[4|). 
The pressure dependence of AC sets constraints on any 
theory of the HO involving quantum criticality. 

In conclusion, using Larmor diffraction to determine 
the lattice constants and their distribution at high pres- 
sures and with very high precision, we were able to show 
that the SMAF in URu2Si2 is purely parasitic. Moreover, 
the HO and LMAF are separated by a line of first-order 
transitions ending in a bicritical point. This is the be- 
havior expected for the HO being, e.g., incommensurate 
orbital currents, helicity order, or multipolar order. 
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